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OBLIQUE-SHOCK RELATIONS AT HYPERSONIC SPEEDS 
FOR AIR IN CHEMICAL EQUILIBRIUM 
By W. E. Moeckel 


SUMMARY 

Oblique -shock relations for air in chemical equilibrium have been 
calculated for flight velocities up to 25,000 feet per second at altitudes 
up to 200,000 feet. Results show that those shock parameters which are 
functions only of Mach number normal to the shock for an ideal gas are 
strongly influenced by flight altitude (initial conditions), as well as 
normal Mach number, when dissociation takes place. 

The variation of flow-deflection angle with shock angle differs sig- 
nificantly from that of an ideal gas. At an altitude of 100,000 feet and 
a flight speed of 25,000 feet per second, for example, the wedge detach- 
ment angle is about 14° larger than that obtained for an ideal gas. 


INTRODUCTION 

Tables and charts of oblique-shock relations for air are currently 
available only for downstream temperatures less than 5000° R (ref. l) . 
Recent computations by the National Bureau of Standards of the properties 
of air in chemical equilibrium at temperatures up to 15,000° K (ref. 2) 
permit extension of these charts to the regions of interest for hypersonic 
aerodynamics of nonslender bodies. Some representative computations, 
based on the data of reference 2, have recently been completed at the NACA 
Lewis Laboratory and are reported herein. Although these computations are 
not as detailed as those of reference 1, they contain sufficient informa- 
tion to calculate with reasonable accuracy the inviscid flow quantities 
downstream of oblique shocks of arbitrary strength for flight speeds up 
to 25,000 feet per second and altitudes up to 200,000 feet. The principal- 
limitation in the use of the oblique-shock relations presented herein is 
the assumption that the air is in equilibrium downstream of the shock. 

Two evaluations of shock relations for air in equilibrium have been 
published recently (refs. 3 and 4). Both references are limited to 
normal-shock relations. Reference 3 was published before thermodynamic 
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data based on the corrected dissociation energy of nitrogen (9.76 electron 
volts) became available. Reference 4 uses the corrected thermodynamic 
data and compares results with the values obtained with the previously ac- 
cepted nitrogen dissociation energy of 7.37 electron volts. 


ANALYSIS 

In the nomenclature shown in the following sketch, the equations of 



continuity, conservation of momentum, conservation of energy, and state 
are as follows: 

p 2 u 2 = (1) 

2 2 

g(p 2 - Pi) = Pl u l - P2 U 2 ( 2a ) 

v 2 = v-^ (2b) 

2gj(h 2 “ h 1 ) = uj - u 2 (3) 


p 2 - 


RJ 


P2 t 2 

m 2 


(4) 


All symbols are defined in the appendix. 

Equations (l) and (2a) yield the following expression for a pressure 
coefficient based on normal Mach number M^ sin 0: 

p 2 

— - 1 

p l = i _ ^2 

YpM 2 sin 2 0 u l 


( 5 ) 
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Equation (3) can be put into the form of an enthalpy coefficient based 
on normal Mach number as follows: 


^2 


r— - 1 


n 


- 1 


= i 


o 9 
MJ sin 0 



( 6 ) 


In equation (6) it is assumed that ambient temperature is low enough that 
the following ideal-gas relations hold: h n = -,t n and y, = . 

1 PA 1 1 

The relation between flow-deflection angle 8 and shock angle 0 
can be expressed as 


u. 


*2 ^2 u 2 u 2 
tan (0 - 8) = — = — = — — = — tan 6 


or 


tan (0 - 8) = ^2 
tan 0 


(7) 


Equations (5), (6) , and (7) show that all quantities downstream of the 
shock (except the speed of sound) can be determined if the ratio of normal- 
velocity components across the shock Ug/u-^ i s known. For, with Pg and 

hg known, the other state quantities Pg, tg, and mg can be obtained 

from the computed properties of air in equilibrium presented in reference 
2. The determination of these state relations is greatly simplified by 
use of thermodynamic charts for air prepared from the data of reference 2. 

The procedure used to solve for ug/u^ was as follows: 

(1) Compute p 2 and hg from equations (5) and (6), respectively, 
for several assumed values of ug/u-^ near the expected correct value. 

(2) With the values of Pg and hg, use thermodynamic charts to ob- 
tain tg and nig. Compute pg from equation (4). 

(3) Plot Pg/pi from step (2) against assumed u-^/ug values. The 
point of intersection of this curve with the line Pg/Pi = u l/ u 2 yields 
the correct value of Ug/u^. 
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Speed of sound downstream of shock . - The speed of sound of a gas 
in equilibrium is given by 


a 


2 



( 8 ) 


where the subscript S indicates that the derivative is taken along an 
isentrope. Alternatively, the speed of sound can be expressed by 


a 2 = Tg£ = rgJRt 
p m 

where r is the local isentropic exponent 1 , which is given by 

v = /c> log p \ 

\ d lo s P/s 


( 9 ) 


( 10 ) 


From the equation of state, the following equivalent formulas are derived: 



The variations of log p and log p with log ^ along isentropes 

were obtained from thermodynamic charts plotted from the data of reference 
2 and from the equation of state. Values of r were then obtained from 
equations (ll) and (12) by measuring the local slopes of these curves. 
Results are shown as a function of enthalpy in figure 1 for several val- 
ues of S/Rq« Although some interpolation error is involved in reading 
the thermodynamic charts, the values of r obtained from equations (ll) 
and (12) agreed within ±0.02. 

With r known, the speed of sound is calculated from equation (9) 
after p and p are determined. 


1 The isentropic exponent is, of course, not equal to the ratio of 
specific heats when dissociation takes place. Definition of r as in 
eq. (10) avoids the necessity of calculating specific heats for the dis- 
sociating gas . 
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RESULTS 

Computations were carried out for the following combinations of ini- 
tial velocity and altitude: 


Altitude, 

Initial velocity, 

ft 

ft /sec 

100,000 

5,000 

10,000 

15.000 

20.000 
25,000 

50,000 \ 

f 5,000 

150,000 J 

< 15,000 
L 25,000 

200,000 

25,000 


The ambient conditions used at each altitude are given in table I. 

The computations are simplified considerably if it is noted that 
u g/ u l is, for given initial conditions, a function only of the normal 
Mach number M^ sin 0. (in general, if the initial enthalpy is suffi- 
ciently large that hp cannot be set equal to Cp^t^, then ug/u^ is 
a function of the initial normal velocity u^ instead of Mj_ sin 0.) 
Consequently, the pressure ratio and enthalpy ratio for given initial 
conditions are also functions only of Mp sin 0 . Since hg and pg 

determine Sg, tg, mg, and pg, these quantities are also functions only 
of Mp sin 0 . 


Normal-Velocity Ratio ^/up 

The ratio ug/up is plotted in figure 2 as a function of the normal 
Mach number Mp sin 0. Shown for comparison is the variation obtained 
from the constant y formula: 


^2 

U 1 


r + i| 



l + 


(13) 
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For normal Mach numbers less than 4, the real-gas results differ 
little from the y = 1.4 values. In the range of normal Mach numbers 
from 4 to 25, the dissociation of air becomes more and more significant, 
and the normal- velocity ratio approaches the r = 1.15 curve. The effect 
of altitude (initial temperature and pressure) is quite appreciable. The 
initial divergence from the y = 1.4 curve is due to oxygen dissociation. 
For Mj_ sin 0 near 14, the oxygen is almost entirely dissociated, and 

the curves tend toward the horizontal. For somewhat larger values of 
sin 0, the dissociation of nitrogen becomes significant, and the 

curves again move downward. 


Flow Deflection Across Shock 

The relation between flow-deflection angle across the shock 6 and 
shock angle 0 is shown in figure 3(a) for an altitude of 100,000 feet. 
The constant y curves, for y = 1.4, seriously underestimate the deflec- 
tion angle for shock angles greater than about 30° when the flight speed 
is greater than 5000 feet per second. The detachment angle, for example, 
is 60.1° for V-l = 25,000 feet per second, as compared with 45.6° for in- 
finite Mach number with y = 1.4. The effect of altitude is shown in 
figure 3(b) . 


Mach Number Behind Shock 

In figure 4, the Mach number downstream of the shock is plotted 
against shock angle for several flight speeds at an altitude of 100,000 
feet. Of particular interest is the fact that the sonic point moves from 
0 = 67° at 5000 feet per second to 0 = 75° at 25,000 feet per second. 
(For y = 1.4, the sonic point is 0 = 68° for M^ = °».) 


Mass -Flow Ratio 

The ratio of the mass flow per unit area across the shock P‘^ 2 / Pi^i 

is shown in figure 5. This ratio is the reciprocal of the stream-tube- 
area ratio across the shock. The maximum mass-flow ratio occurs at shock 
angles considerably smaller than that at the sonic point. This result is 
due to the strong entropy change with shock angle and is also obtained 
using ideal-gas theory. 


Temperature 

The temperature downstream of the shock is shown in figure 6 as a 
function of normal Mach number. Temperatures are higher at 50,000 feet 
than at 100,000 feet because of the higher pressure at lower altitudes 
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which inhibits dissociation. The initially, higher downstream tempera- 
tures at 150,000 feet are, of course, due to the higher ambient tempera- 
ture at this altitude. 


Entropy Increase 

The increase in entropy across oblique- shock waves is shown in fig- 
ure 7 as a function of normal Mach number. Considerably larger increases 
in entropy are obtained at high normal Mach numbers for the real gas than 
for the ideal gas. 


Other Parameters 

With the oblique-shock relations plotted in figures 2 to 7, all 
other flow parameters downstream of the shock are easily obtained. The 
density ratios across the shock are simply the reciprocals of the normal- 
velocity ratios of figure 2. Pressure and enthalpy downstream of the 
shock are also obtainable from figure 2 with the aid of equations (5) and 
(6). The absolute velocity and speed of sound downstream of the shock 
can be calculated from 


or 


and 


•w-2 — 2 i tt 2 

v 2 - u 2 + V 1 


cos 2 0 



(14a) 


(14b) 


a 2 " M 2 

where u 2 / u i is obtained from figure 2, 2/(^1 

M 2 from figure 4. 


(15) 


from figure 5, and 
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APPENDIX - SYMBOLS 

speed of sound, ft/sec 

specific heat at constant pressure 

specific heat at constant volume 

acceleration due to gravity, 32.16 ft/sec 2 

enthalpy, Btu/lh (h * 0 for undissociated air at t = 0) 

mechanical equivalent of heat, 778 ft-lb/Btu 

Mach number 

molecular weight, lb/lb-mole 

molecular weight of undissociated air (28.86 for composition of ref. 

2 ) 

pressure , Ib/sq ft 

gas constant, 1.987 Btu/(lb-mole) (°R) 

R/mo 

entropy , Bt u/ ( lb ) ( °R ) 
temperature, °R 

reference temperature of reference 2, 491.69° R 
velocity component normal to shock, ft/sec 
resultant velocity, ft/sec 

velocity component parallel to shock, ft/sec 
isentropic exponent (y-^ = 1.4) 
flow-deflection angle across shock, deg 
shock angle, deg 
density, lb/cu ft 
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Subscripts : 

S differentiation along isentrope 

1 conditions upstream of shock 

2 conditions downstream of shock 
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TABLE I. - ASSUMED AMBIENT CONDITIONS 


Altitude, 

ft 

Pressure, 

*i' 

atm 

Density, 

p i> 

Tb/sq. ft 

Temperature , 
° R 

Dimensionless 
enthalpy , 

h l/ R O t O 

Dimensionless 

entropy, 

Si / R 0 

0 

1.0 

0.0765 

519 

3.70 

23.6 

50,000 

0.1145 

0.01165 

392,4 

2.795 

25.0 

100,000 

0.0106 

0.001065 

392.4 

2.795 

27.4 

150,000 

0.00142 

9 .81 x 10 -5 

573.5 

4.085 

30.9 

200,000 

0.000314 

1.988X10' 5 

619.4 

4.41 

32.7 


Isentropic exponent 




1.40 


1.32 


1.24 


1.16 


1.08 
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Figure 1. - Isentropic exponent for air In equilibrium. 
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Normal-velocity ratio, u 2 /’ 
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Figure 2. - Normal-velocity ratio across shock waves as function of normal Mach number 
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(a) Effect of velocity. Altitude, 100,000 feet. 
Figure 3. - Shock angle as function of flow-deflection angle. 
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(b) Effect of velocity and altitude. 

Figure 3. - Concluded. Shock angle as function of flow -deflect! on angle 
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Mach number behind shock 
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Shock angle, 0 , deg 

Figure 4. - Mach number downstream of shock. 
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Figure 7. - Entropy increase across shock waves. 
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